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The magnon dispersion in Ca 2 Ru 04 has been determined by inelastic neutron scattering on 
single crytals containing 1% of Ti. The dispersion is well described by a conventional Heisenberg 
model suggesting a local moment model with nearest neighbor interaction of J=8meV. Nearest 
and next-nearest neighbor interaction as well as inter-layer coupling parameters are required to 
properly describe the entire dispersion. Spin-orbit coupling induces a very large anisotropy gap 
in the magnetic excitations in apparent contrast with a simple planar magnetic model. Orbital 
ordering breaking tetragonal symmetry, and strong spin-orbit coupling can thus be identified as 
important factors in this system. 


The properties of strongly correlated systems with sig¬ 
nificant spin-orbit coupling (SOC) present a challenging 
problem. The intensively studied example is the reduc¬ 
tion of the magnetic state of Ir^+ (electronic structure 
or t 2 g, Leff = 1, S=l/2) to an effective Kramers 
doublet with J=l/2 [l|. But even more drastic effects 
can be expected for heavy ions with d'^ occupation 
Leff=l, S=l), e.g. in Ir®+, Ru'*+, Os‘*+ etc.Q. Accord¬ 
ing to Hand’s rules (generalized for ions sensing crystal 
electric fields) the ground state should be a nonmagnetic 
singlet with J=0, see e.g. reference ii- And indeed iso¬ 
lated Ir^+ ions and also most of concentrated Ir®+ com¬ 
pounds are nonmagnetic, although a few magnetic Ir®+ 
cases are known Q. In a solid magnetic order can occur 
even if the ground state of an isolated ion is a singlet, 
see chapter 5.5 in reference 0], but it requires a strong 
exchange interaction, so that the exchange splitting of 
excited magnetic states (in the Ru'^+ case a J=1 triplet) 
is larger than the energy difference between the ground- 
state singlet and the excited triplet, which is given by 
the SOC parameter A. The SOC can also be at least 
partially suppressed by a non-cubic crystal field (CF), 
Anoncub, which splits the Og {Lgff=l) triplet and stabi¬ 
lizes real orbitals. Both these factors, CF and magnetic 
interaction, can combine to suppress the J=0 state and to 
eventually induce the magnetically ordered ground state. 
In terms of energy scales, one should expect such mag¬ 
netic ordering for Anoncub + R ■ Hexch > A, which seems 
quite unlikely for Ir^+ where A=^ = C/2 amounts to 0.2 
to 0.25 eV (C is the atomic spin-orbit parameter). But 
for 4(i compounds this relation can easily be reached, as 
for Ru^+ A ~ 0.075 eV[3,|^. Indeed, practically all Ru"‘+ 
compounds order magnetically aside from the metallic 
ones - and even some metallic ruthenates are magnetic, 
such as the ferromagnetic metal SrRuOa. The persisting 
role of SOC in these magnetic Ru^"*" compounds is an 
intriguing open issue. 

Ca 2 Ru 04 (CRO) is such a Ru^+ case, which has 
been intensively studied as the Mott-insulating ana- 
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FIG. 1: Intensity distribution in energy versus scattering vec¬ 
tor, Q, planes taken at 2K around the (1.5 0.5 0) magnetic 
zone center, a) and c) show the symmetrically equivalent dis¬ 
persion along (0,^,0) and (^,0,0) direction, subplot b) and d) 
along (^,^,0) and (^,-^,0). The color coding corresponds to 
the raw data. Open symbols indicate the dispersion obtained 
by fitting single scans. Data were taken with final energies 
of 35meV for constant Q scans at high energy transfer and 
14.7 meV elsewhere. 


logue of the unconventional superconductor Sr 2 RuO4 0 - 
[l^. CRO exhibits a metal-insulator (MI) transition at 
357 K which is accompanied by a flattening of the RuOe 
octahedrafil-fl^ . This flattening continues upon further 
cooling until it saturates near the onset of magnetic or¬ 
der at Tat^IIOK. The magnetic structure is antiferro¬ 
magnetic (AFM) with moments aligned parallel to the 
layers [T2|- The electronic structure has been stud- 
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ied by various approaches |13l - ll7| . From the spectro¬ 
scopic study of CRO it was concluded that SOC indeed 
plays an important role but is not sufficiently strong to 
stabilize the J=0 state [l^. Density functional theory 
calculations indicate a pronounced shift in orbital polar¬ 
ization leading to almost full occupation of the dxy levels 
at low temperature [l3 - [l8j| . More recently the J=0 state 
was explicitly proposed for CRO ii. Starting from 
the scenario of strong SOC and including noncubic CF 
and intersite exchange, the magnetically ordered state in 
CRO is reproduced and several unusual features of the 
magnetic excitation spectrum of CRO are predicted, such 
as a peculiar shape and large width. The alternative, 
more conventional picture is to attribute the magnetism 
of CRO to the conventional S^l state of Ru^"*" ions, with 
SOC playing a less significant but still prominent role. In 
this case one can describe the magnetic state, including 
spin waves, by the usual exchange Hamiltonian. 

Here we present an inelastic neutron scattering (INS) 
study and spin-wave calculations of the magnetic excita¬ 
tions in CRO. We find that a conventional model can well 
describe the obtained dispersion, while there are consid¬ 
erable differences with the proposed J=0 modelQ. Most 
interestingly there is a sizable spin-gap which indicates 
that rotating the magnetic moment within the layers 
costs large energy. The breaking of the local tetragonal 
symmetry and the associated orbit al p olarization, which 


has been neglected in theory so far |14l - ll8| , are important 


parameters to understand the magnetism in CRO. 

CRO single crystals of of several lOOmm^ volume con¬ 
taining 1% of Ti were obtained by the travelling solvent 
floating zone method in a mirror furnace (Canon SCI- 
MDH11020-CE). We added 1% of Ti as this seems to 
avoid bursting of the crystal upon cooling below the MI 
transition. Characterization by magnetic susceptibility 
and by neutron diffraction experiments indicate a mag¬ 
netic transition at Tjv=112 K and no significant impact 
of the very small amount of Ti. INS experiments were 
performed with the INS triple-axis spectrometer. Most 
experiments were performed with fixed final momentum 
of A:/=2.662 A-i (E/=I4.7 meV); some scans at high en¬ 
ergy transfer or aiming at better resolution were per¬ 
formed with kf=i.l and A:/=I.97A“^, respectively. We 
studied the magnon dispersion in the two scattering ge¬ 
ometries (100)/(010) and (110)/(001) in reduced units 
of the tetragonal lattice [l2|. Eor both setups two crys¬ 
tals were coaligned. The sample of the second scattering 
plane was essentially untwinned as determined on the 
INS spectrometer. 

Fig. I shows color mappings of the measured inten¬ 
sity distribution. Due to the weakness of scattering in 
CRO (small moment and rapidly decreasing form factor) 
contaminations by various phonon branches are highly 
visible. By analyzing and comparing results taken in dif¬ 
ferent Brillouin zones and geometries the dispersion can 
be unambiguously determined. Magnon excitations start 
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FIG. 2: Several characteristic scans: a) Constant energy scans 
at (1.5,fc,0) taken at 2 K fitted with gaussians and background, 
b) Phonon scans taken at Q=(^,2,0); the lines correspond to 
the folding of the resolution function with a simple linear 
phonon dispersion. No additional parameter is needed to de¬ 
scribe the shape of the intensity profile, c) and d) show the 
energy scans at the zone center (1.5 0.5 0) and at (0.5 0.5 1). 


at the AFM Bragg points 2 n|+i integer 

n^, n/; and ni. However, there is a sizeable spin gap of 
13.04(5) meV. For a square planar antiferromagnet the 
magnon dispersion extends from Q=(0.5 0.5) to (0.75 
0.75) in [1 1] direction, as (1 1) is a Bragg point, and to 
(0 0.5) in [1 0] direction. Q=(0.25 0.25) and (0 0.5) are 
AFM Brillouin zone boundaries. In CRO there is, how¬ 
ever, a severe structural distortion 0. Some character¬ 
istic scans performed to determine the magnon dispersion 
in CRO are shown in Fig. 2. Constant energy scans at in¬ 
termediate energy cut through the magnon cones at two 
positions. Constant Q scans taken just at the AFM zone 
center show a characteristic asymmetric shape, see Fig. 
2 (c) and (d): Intensity rapidly increases when crossing 
the spin gap and slowly diminishes with further energy 
increase. We have calculated the folding of the spin- 
wave dispersion including its expected signal strength 
with the experimental resolution using the RESLIB [20| 
package and verified that scans across transversal acous¬ 
tic phonons are well reproduced, see Eig. 2 (b). The steep 
spin wave dispersion perfectly describes the asymmetric 
shape of the spectra taken at the zone center, see Fig. 2 
(c). The total width of the dispersion is low, as maxi¬ 
mum energies of 37.8(3) and 41.2(5) meV are reached at 
the magnetic zone boundaries, (0.5 0 0) and (0.25 0.25 0). 

In order to describe the magnon dispersion we use 
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FIG. 3: (Dispersion of the magnon branch along the main 
symmetry directions at T=2 K. The open symbols indicate 
the values obtained by fitting the raw data scans with Gaus- 
sians or by folding the resolution function with the modeled 
dispersion. Lines correspond to the spin-wave calculations 
with the Heisenberg model as described in the text. 


a conventional Heisenberg model with a single-site 
anisotropy term arising from SOC: H = j JijSj • Sj — 
We set S=0.67 following the neutron diffrac¬ 
tion study [^. The sum runs over pairs of magnetic ions, 
so that each pair or bond appears twice. Spin waves 
were calculated with the Holstein-Primakoff transforma¬ 
tion as described in references [U . We include the 
nearest-neighbor magnetic exchange of J=8meV, next- 
nearest neighbor interaction along the orthorhombic a 
and b directions of Jn„a=Jrmb=0.7meV and an AFM 
coupling between neighboring layers. The next-nearest 
neighbor interaction is chosen isotropic, as the twinned 
crystal used in the (1 0 0)/(0 1 0) geometry prohibits dis¬ 
tinguishing these directions. The need for the additional 
parameter can be seen when comparing the magnon en¬ 
ergies at q=(0.25 0.25 0) and (0.5 0 0), which are identi¬ 
cal in the model with only nearest-neighbor interaction. 
The interlayer coupling, Jc=0.03meV is the only param¬ 
eter that breaks the tetragonal symmetry in our model 
aside from the single-ion anisotropy. Note, however, that 
the crystal structure is orthorhombic lifting the degen¬ 
eracy of magnetic interaction parameters. We chose the 
AFM interaction between the Ru at (0,0,0) and that at 
(0,0.5,0.5) (in the orthorhombic cell [9|), which stabilizes 
an A centered magnetic structure with magnetic space 
group Pbca 0. 

The magnetic moment in CRO points along the or¬ 
thorhombic b direction. Therefore, one might expect a 
large gap for the magnetic excitations involving rota¬ 
tions of the moment out of the Ru 02 layers, and much 
softer in-plane modes. The latter are described by the 
expectedl y s mall in-plane anisotropy. Following refer¬ 
ence [13, I 23 II both branches can be described simulta¬ 
neously with two anisotropy parameters. Surprisingly, in 
CRO the in-plane anisotropy turned out to be extremely 
strong. The magnon dispersion starts at 13.04(5) meV 
which we may identify with the in-plane gap. There is no 
magnon branch at lower energy as is clearly shown in the 


intensity maps. We find some weak intensity at an energy 
transfer of 5meV appearing near Q=(1.5 0.5 0), but this 
signal is flat in energy and restricted to the AFM zone- 
center. Furthermore, this signal is much weaker than 
the magnon modes at higher energies suggesting a pos¬ 
sible origin associated with the Ti impurities and some 
local effect. As shown in Fig. 2 (d) and 3 there is a 
finite inter-layer dispersion visible in the scans taken at 
Q=(0.5 0.5g/) with the second untwinned crystal. The 
tetragonal [110] direction corresponds to orthorhombic 
b in the used mounting and thus to the direction of the 
magnetic moment; therefore, the transverse magnon with 
in-plane polarization (thus parallel to orthorhombic a) 
fully contributes. Also in the other configuration there is 
a clear difference in spectra taken at Q=(0.5 0.5 0) and 
=(1.5 0.5 0). For the twinned sample we superpose AFM 
zone-centers and zone boundaries. For a twinned crys¬ 
tal c polarized magnons will always contribute, while for 
the in-plane magnon the geometry condition that only 
magnetic components perpendicular to Q contribute, 
suppresses some modes. The fact that we see a clear 
difference at various ( 2 n/^+i unambiguously 
shows that the modes dispersing between 13.04(5) and 
14.2(1) meV posses an in-plane polarization. This fur¬ 
thermore agrees with the Qi dependence of the signal. 
We may thus conclude that the lowest magnon branch in 
CRO possesses an in-plane character and that it disperses 
between 13.04(5) and 14.2(1) meV along the c direction 
and up to 41.2 and 37.8 meV along (^ ^ 0) and (^ 0 0) 
paths, respectively. We cannot identify the c polarized 
modes as they may remain hidden in the asymmetric 
shape. There is some evidence for a nearly flat branch 
around 36meV, but we cannot fully rule out that these 
modes are purely nuclear or that they carry longitudinal 
polarization. For simplicity, the experimental dispersion 
is described by an easy-axis anisotropy [ 2 ^, see Fig. 3. 

The magnon dispersion including its large gap can be 
very well described within the spin-wave theory suggest¬ 
ing a conventional local moment S'^l magnetism with a 
strong - but not decisive - impact of SOC. Starting from 
the other scenario, a spin-orbit driven J=0 singlet nature 
which is rendered magnetic by noncubic CF and intersite 
exchange, Akbari and Khaliullin predicted several un¬ 
usual features of the magnetic excitation spectrum, such 
as the energy continuously softening from the value A 
at F, and the presence of extra modes in some part of 
the spectrum. Our results, however, do not support this 
model 0. First, the observed dispersion is much flatter 
than this prediction, as it does not reach energies of the 
order of the expectedly large values of A, and as there is a 
strong gap. Second, the singlet picture predicts a contin¬ 
uously increasing dispersion near the AFM zone bound¬ 
aries, while our experiments find the saturation predicted 
by the Heisenberg model, see Fig. 1 and 2. The Heisen¬ 
berg scenario also implies several branches: two transver¬ 
sal branches arise from the orthorhombic anisotropy (in- 
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FIG. 4: Temperature dependence of the magnetic scatter¬ 
ing at the AFM zone center (1.5,0.5,0) measured with the 
Si monochromator and fcF=1.97 (a) and 2.662A“^ (b). The 
sharp peak associated with the in-plane spin gap mode softens 
and broadens considerably. The inset shows the energy of the 
spin gap scaled to its low temperature value (blue) compared 
to that of a Raman signal taken from (red). 


plane and c polarized), and longitudinal modes can exist 
in CRO being on the border to itinerancy. 

Using the standard description, with the hopping pa¬ 
rameters t ^100 meV, obtained by ab-initio calculations 
1^ 241, and using the Hubbard’s U'^2eV, we would ob¬ 
tain for the exchange constant J=2t^/U ^lOmeV, in 
good agreement with our experimental finding. How¬ 
ever, CRO is not a strong Mott insulator with completely 
localized electrons as it is already indicated by the low- 
lying MI transition. In this case the basic J=0 Ansatz 
may be not a good starting point, as the J=0 state can 
be suppressed by electron hopping. Also for Ir^+ (specifi¬ 
cally for Na 2 lr 03 ) the sizable hopping modifies the whole 
picture [2^, [2^ , leading to novel quasi-molecular orbital 
states with reduced impact of SOC. The conspicuous but 
typical absence of J=0 physics in most of the Ru‘^+ ma¬ 
terials seems largely connected with the hopping. 

Another argument in favor of the applicability of the 
usual picture of Ru^"'' ions (S^l) is the strong flatten¬ 
ing of RuOe octahedra [llj occurring below the MI 
transition. Such distortion is typical for the usual Jahn- 
Teller effect: it stabilizes the doubly-occupied dxy orbital, 
leaving two electrons on dxz and dyz- In such state the 
orbital moment and spin-orbit interaction are partially 
quenched. The sign of this distortion proves that in this 
system the Jahn-Teller effect is stronger than the SOC 
which would have caused the opposite distortion and CF 
splitting 0. Recent spectroscopy data [13 confirm this 
significant splitting of t 2 g orbitals. 

On the other hand the observation of the strong in¬ 
plane magnetic gap is remarkable for a layered system. 
It underlines the relevance of the SOC in CRO even in 
the conventional scenario. Several Raman scattering ex¬ 
periments observed an additional signal in Big symmetry 
appearing in the AFM phase 27H29|. This feature was 
interpreted as a two-magnon excitation, but our results 


clearly rule out such explanation. The Raman feature 
appears at 102cm“^=12.6meV at 10 K which is much 
below the energies for two magnon excitations and the 
expected peak in the two-magnon density of states (near 
SOmeV). Instead this energy agrees with that of the in¬ 
plane gap mode. The single magnon mode, however, is 
not Raman active in first approximation, but SOC can 
induce a finite signal. The temperature dependence and 
the extreme broadening of the Raman signal at higher 
temperature agree reasonably well with the correspond¬ 
ing behavior of the magnon gap, see inset in Fig. 4 (b). 

The magnetic in-plane anisotropy in CRO must origi¬ 
nate from SOC and from an orbital arrangement break¬ 
ing tetragonal symmetry. There have been many exper¬ 
imental and theoretical analyzes l3“18| elucidating the 
change of the orbital polarization upon cooling and the 
increasing occupation of the dxy versus the dxz/dyz or¬ 
bitals following the flattening of the RuOe octahedron. 
This distortion possesses Eg symmetry, which is the most 
frequently analyzed in Jahn-Teller models [s^- The t 2 g 
orbitals, however, also couple to the T2g octahedron dis¬ 
tortions (33 which break tetragonal symmetry in the 
case of CRO but which were neglected so far. The tem¬ 
perature dependence of the crystal structure of CRO in 
the insulating phase reveals an ongoing elongation of the 
RuOe octahedra iEl along the orthorhombic b direc¬ 
tion along which moments align. This distortion corre¬ 
sponds to the T 2 g ’’scissor” mode of the free octahedron 
[SOj lifting the dxz/dyz degeneracy. Similar to a tetrag¬ 
onal distortion, e.g. along z-axis, which would stabilize 
orbitals with E = ±1, {dxz ± idyz), and which, by SOC 
would orient spins along z-direction, (or a trigonal elon¬ 
gation along [111] (in cubic setting), which would make 
[111] an easy axis, see e.g. 0), such T2g distortion (elon¬ 
gation along b axis) makes the orthorhombic b-direction 
the easy axis. 

In conclusion we have studied the magnon dispersion in 
CRO which considerably differs from recent predictions 
for a J=0 singlet ground state. Instead, the dispersion 
is well described in a local moment Heisenberg model 
with strong anisotropy terms yielding a nearest-neighbor 
exchange interaction of J=8meV which agrees with the 
large calculated hopping integrals. Large hopping seems 
to be the main cause for the suppression of the J=0 state 
in Ru^“'" compounds. On the other hand, the remarkably 
strong in-plane anisotropy clearly shows that considering 
tetragonal crystal fields is insufficient. There is impor¬ 
tant orbital polarization breaking tetragonal symmetry, 
which is related to the prominent elongation of RuOg 
octahedra along the orthorhombic b direction and which 
renders spin-orbit coupling still active in this system. 
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